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Abstract. Based upon the splendid experimental result of Gregor et al. ]S\, 
we present a dynamical system model on the motion of individual amoebae. 
The dynamics of our model qualitatively agrees well with Gregor et al.'s ex- 
f^ ^ perimental result on the development of collective behavior in social amoebae. 
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r^ . 1. Introduction 

O^. Complexity is the current novelty in Science. Chaos as a part of complexity 

has proved to be an important concept. Besides chaos, another concept that is very 
elusive in current complexity pursuits is emerging '5^ . Among the phenomena that 
^ ' may be termed emerging, collective behavior of individual elements is the most in- 

^^ . triguing. The potential importance of collective behavior was stressed long ago [1]. 

Recent more advanced experiments on microorganisms have revealed impressive 

^^ i collective behavior via dynamical quorum sensing j2j. We are especially impressed 

by the elegant experiment of Gregor et al. [3] , which can manipulate and monitor 
an individual amoeba in the social amoebae Dictyostelium discoideum. An amoeba 
can synthesize and release into the extracellular space, a chemoattractant named 
cAMP (cyclic adenosine 3',5'-monophosphate). When it is stimulated with extra- 
cellular cAMP, an amoeba responds by synthesizing and releasing more cAMP [7] . 
The extracellular cAMP guides the amoebae aggregation to form fruiting bodies. 
The main discovery of Gregor et al. |3] is that collective behavior originates from 
the specific amoeba which is stimulated with extracellular cAMP. Here we shall 
introduce a dynamical system model to simulate the motion of the amoebae aggre- 
gation. Using mathematical equations to describe biological or chemical phenomena 
has a long history [4] [5] [6] , and has established itself an important scientific area. 

2. The model 

Now we shall describe our model. Let N be the total number of amoebae 
{N = 180 in [3"). Denote by Xn the location of the n-th amoeba in the space. 
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The space can be S-dimensional, 2-dimensional, or l-dimensional. The dynamical 
system modehng the motion of the amoebae is given by 

(2.1) itn = -^ Va;„G(|i;„ -i„J) - ^ Vx„H{\Xn -x„J) + Snn*F{t), 

jeJ keK 

where '■' denotes time derivative d/dt, \ \ denotes the Euchdean vector length, 
V denotes the gradient, G is an attractive potential, H is an attractive-repulsive 
potential, J and K represent the neighbors, the G and H terms represent the 
forces of the extracellular cAMP, F{t) is the stimulating force exercised by injecting 
extracellular cAMP at a specific amoeba, and 

5nn* — 1 for n = 11*, =0 for n ^ n*] 

where n* denotes the specific amoeba. We are interested in the following examples 
for G, H and F: 

Oi „ 



G{r) = 
H{r) - 
F{t) = 
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where (a,/3,7,a, 6,c, cr) are positive constants. The model now takes the specific 
form 
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(2.2) 



cosi. 



The setup for the numerical simulation of (|2.2p is as follows: The space is the 
2-dimensional plane, 

a = 50, a = 3, /3 = 1, 6 == 2, 
c=l, cr^O.l, 7 = 2, A^=100; 
thus each Xn is a 2- vector. For any n, 

if \xn - i„ J < 4, then j e J; 
similarly 

if \xn — Xrik I f^ 4, then k E K. 

The initial locations of the 100 amoebae are uniformly distributed in a circular 
chamber of radius 10. The initial velocities of the 100 amoebae are uniformly 
distributed in 

[-0.001,0.001] X [-0.001,0.001]. 

We use the Runge-Kutta scheme with time step length 0.01 to conduct the nu- 
merical simulation. The motion of the 100 amoebae modeled by (|2.2p qualitatively 
agrees well with the experimental result of Gregor et al. [3], see Figures [T]|ni 
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Figure 1. The locations of the amoebae at time t = 0. 
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Figure 2. The locations of the amoebae at time t = 100. 
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Figure 3. The locations of the amoebae at time t = 200. 
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Figure 4. The locations of the amoebae at time t = 300. 



COLLECTIVE BEHAVIOR IN SOCIAL AMOEBAE 




-15 -10 -5 5 10 15 



Figure 5. The locations of the amoebae at time t = 400. 
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Figure 6. The locations of the amoebae at time t — 500. 



